Exploring the order parameter symmetry of p-wave Fermi condensates 
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We discuss the time-of-flight expansion of dilute p-wave Fermi condensates on the BEC side of 
Feshbach resonances, as a way to extract information about the order parameter symmetry for su- 
perfluidity. We show that the cloud profile is in general sensitive to the interaction strength between 
fermions, the magnitude and direction of external magnetic fields, and to the angular momentum 
projection of the order parameter. In particular, due to the anisotropic nature of p-wave interactions 
we show that the time-of-flight expansion of a p-wave superfluid is anisotropic even if the superfluid 
is confined to a completely isotropic trap, unlike the case of Bose or s-wave Fermi condensates, which 
under the same circumstances expand isotropically. Furthermore, we demonstrate that expanding 
p-wave superfluids released from axially symmetric traps experience anisotropy inversions, where 
the aspect ratio between the axial and radial directions change during expansion, while the axial 
symmetry can also be lost reflecting the spatial anisotropy of the underlying interaction. 

PACS numbers: 03.75.Ss, 03.75.-b, 05.30.Fk 



One of the next frontiers of research in ultra-cold atoms 
is the search for superfluidity in higher angular momen- 
tum channels, which are quite rare in nature. With the 
exception of liquid 3 He, there are no known neutral su- 
perfluids with angular momentum in the p-wave channel 
or higher. In addition, there are very few known ex- 
amples of p-wave and d-wave superconductors (charged 
superfluids) in standard condensed matter physics, where 
the degree of control of the interactions is not existent. 
In contrast, for ultra-cold atoms, interactions in higher 
angular momentum channels can be tuned via the appro- 
priate choice of Feshbach resonances. 

In general, a key point for the understanding of su- 
perfluidity is the determination of the order parameter 
symmetry, and some ingenuity is usually necessary to 
extract this information experimentally in standard con- 
densed matter systems. However, for ultra-cold atoms, 
time-of-flight measurements can provide valuable infor- 
mation about their superfluidity and the corresponding 
order parameter symmetry, as discussed below. 

Timc-of-flight experiments are powerful methods to 
study ultra-cold atoms, and have been applied to confirm 
the existence of Bose-Einstein condensation (BEC) and 
to explore superfluid dynamics of Bose systems [l|, |2|, 13, \4 , 
[B|, [||. Furthermore, this type of measurements were used 
to probe signatures for superfluidity of fermions with s- 
wave interactions 0, 0, @ . 

Timc-of-flight techniques were also used to detect p- 
wave Feshbach resonances of ultra-cold fermions in har- 
monic traps @, H, 0, E3] and in optical lattices 11 1, 
and paved the way for the exploration of p-wave Fermi 
superfluidity, where fermions are paired in the angular 
momentum channel 1=1. Using these observed reso- 
nances ,j>-wave molecules have been produced in the labo- 
ratory 111, EES m 

for fermion isotopes 40 K and 6 Li, bring- 



ing the experimental community a step closer to realizing 
p-wave Fermi superfluidity of ultra-cold atoms. In addi- 
tion, other Fermi atoms like 173 Yb have also been cooled 
to quantum degeneracy [15[ , and are likely to possess p- 
wave Feshbach resonances, thus permitting the study of 
p-wave molecules or superfluids in non- alkali fermions. 

In contrast to their s-wave counterparts, p-wave Fes- 
hbach resonances have characteristic features including 
splitting of peaks depending on hyperfine (pseudospin) 
states as in 6 Li and 40 K (i.e., |11), |12) + |21), and 
|22)) @, i, @, OH [HI, and splitting of peaks depend- 
ing on angular momentum projections as in 40 K (i.e ., 
the magnetic quantum number mg = or ±1) (jl. [llj]. 
Hence, these splitting may allow the separate tuning 
of p-wave scattering parameters in different pseudospin 
and/or mg channels, such that p-wave interactions may 
be anisotropic in both pseudospin and angular momen- 
tum states. 

Theoretical studies of p-wave systems focusing on the 
evolution from Bardeen-Cooper-Schrieffer (BCS) to BEC 
superfluidity have emphasized the existence of quantum 
phase transitions However, these works dealt 

only with thermodynamic properties, while here we study 
the time-of-flight expansion of ultra-cold fermions with 
p-wave interactions and focus on the BEC side of Fes- 
hbach resonances, where the critical temperatures are 
higher. Our main conclusions are as follows. We show 
that the cloud profile reveals the anisotropic nature of 
the p-wave interaction and provide valuable information 
about the order parameter for p-wave superfluidity. In 
particular, we find that a completely isotropic cloud ex- 
periences anisotropic expansion driven by p-wave inter- 
actions, in contrast to the isotropic expansion of Bose 
and s-wave Fermi systems under the same circumstances. 
For a cigar-shaped cloud, we find that the aspect ratio 



2 



between the axial and radial directions of a p-wave Fermi 
condensate undergoes an anisotropy inversion driven by 
the anisotropy of the harmonic potential, like in Bose 
and s-wave Fermi systems. However, the axial symmetry 
of p-wave Fermi condensates may be broken during the 
expansion reflecting the spatial anisotropy of the under- 
lying interaction, in sharp contrast to Bose and s-wave 
Fermi systems, where the axial symmetry is always pre- 
served. 

To describe the time-of-flight expansion of p-wave 
Fermi condensates, we consider dilute fermions with mass 
to in a single hyperfine state (pseudospin) . The time- 
dependent Hamiltonian is (with h = Hb = 1) 



H(t) = \ drU + (r,f) 



-± + U ext (r,t) 
2m 



diV(r, t)^ (r', t)V(r - r')^(r', t)V(r, t) 



where ip^ (ip) are creation (annihilation) operators, 
?7 ex t(r, t) is the time dependent trapping potential, and 
V(r — r') denotes the fermion-fermion interaction. The 
superfluid state in a harmonic trap is characterized by 
two length scales: the coherence length £ and the length 
scale Ru over which U cx t varies. In the limiting case 
where the trapping potential varies slowly in compari- 
son to the coherence length (Ru ^> £), a semiclassical 
approximation is applicable such that one can first con- 
sider a dilute Fermi gas in free space and add the har- 
monic potential later. With this approximation, the free 
space Hamiltonian takes the form 

H = J2 C(k)4^k + \ £ V(k, k')£t, q i?k >q , (2 ) 

k k,k\q 

where £(k) = k 2 /(2m) — fi is the fermion disper- 
sion shifted by the chemical potential /i, and = 

^k+q/?^-k+q/2 * s ^he crea ti° n operator of a pair 
with center of mass momentum q and relative mo- 
mentum 2k. Here, the momentum space interaction 
V(lt, k') is the Fourier transform of V(r — r'), and 
can be approximated by a separable form V(k, k') = 
-4TE m ,5 m ,r(fc)r(fc')^i: mf (k)r 1 ^(k'), where T(k) = 
kko/(k 2 + k 2 ) is the symmetry function and Yx >mi are 
spherical harmonics. The momentum scale ko is deter- 
mined by the effective range Rq « 1/fco of the interac- 
tion and the diluteness condition (nR^e 3 <C 1) requires 
(/co/^f) 3 3> 1) where n is the particle density and kp is 
the Fermi momentum. The interaction can be related to 
the scattering volume a p through the T-matrix [Hi ]. 
We introduce the bosonic field operator b jf me = 

Ek r ( fc ) Y i*m f ( k ) , /'i +q/2 ^ k+q/2 , its corresponding aux- 
iliary fiend A* and integrate out the fermions to derive 
the Ginzburg-Landau (GL) effective action 



Sgl = So 



/ j m£,m' £ 
q,f ,m f g 



:(q, f)A* (q, ^)A m ,(q, v)+S A , 



around the trivial saddle point A me = near the transi- 
tion temperature T c , where S4 denotes the quartic terms. 
The long wavelength behavior of the system can be de- 
scribed by the expansion of the static part of C^ 1 in 
powers of q, leading to 



ZT 1 ,(q,0)=a 



m£,,m' t 



^1 + .... (3) 



Similarly, the low frequency time evolution of the system 
can be described by the frequency dependent part of £ _1 
in powers of i/q after analytic continuation iv — » vq + i0 + . 
For simplicity, we choose to investigate the cases where 
the order parameter A is nonzero only for the channel 
mi = or me = ±1, but not in a mixed state of me = 
and rug = ±1. In this case, the tensor cjjf" = c^^ij 
becomes diagonal but is in general anisotropic. 

Returning to space-time coordinates and adding back 
the trapping potential, we obtain the equation of motion 
for a p-wave Fermi condensate 



i—A mi (T,t) 



E 



C* V 2 (7 

me % . u m e 

drri£ 2to d m p 



+2U ext (r,t) + -^\A me (r,t)\ 



A m< (r,t), (4) 



where the coefficients a mi , b me , c me , and d me can be ob- 
tained accordingly [18]. The trivial expansion of the 
normal component is included in the evolution of 5*o, 
and is not discussed here. In the BEC limit, the ra- 
tio c % mi /d mi = 1/2, such that Eq. ([4]) reduces the con- 
ventional Gross-Pitaevskii (GP) form for a dilute gas of 
paired fermions with mass 2to. For simplicity and defi- 
niteness, we consider the case of a static harmonic trap 
U exb (r,t) = Ylj=x,-y,z muj j r j/^ for t < which is turned 
off at t = 0. For t < 0, the system is described by 



MoA m ,(r) 



C X7 



a,. 



^— •* d mt 2m d 



-2C/ ex t(R) 



^|A m ,(r)| 2 



A mt (r), (5) 



where [1q is the effective chemical potential. For 
dominant effective boson interactions the Thomas- 
Fermi approximation leads to \A me (r, 0)| = [d me ((Xo — 
2C/ cxt (r))/6 m J 1 /2 for /iQ > 2C/ oxt (r), and |A mf (r,0)| = 
otherwise. When this approximation fails the initial con- 
dition for the time-of-flight expansion can be obtained by 
solving the static equation (O numerically. Furthermore, 
since the trapping potential is separable in the spatial 

variables rj, a scale transformation Rj 



\J cin e /d me can be introduced to remove the 
anisotropy in the gradient term of Eq. ^ , leading to an 
equation with the conventional GP form. 



and Qj — > utj 
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FIG. 1: Aspect ratio r xz = L x /L z as a function of dimension- 
less time r = uj±t during time-of-fiight expansion of a p-wave 
Fermi condensate in the strong attraction (BEC) limit, for ax- 
ially symmetric cigar-shaped traps with e = 0.1 and e = 0.3. 



For t > 0, we introduce the transformation Rj(t) = 
bj(t)Rj(0), where the scaling factors bj(t) satisfy 



d 2 b 3 (t)/dt 2 = n 2 /[A(t)b 3 (t)}, with A(tl 
and initial conditions &j(0) = 1 1^, 20 ] 



b x {t)b y {t)b z {t) 
Under the ap- 
proximation of collisionlcss hydrodynamics, the time evo- 
lution of the condensate wavcfunction 



A mf (R, t) 



exp[iS(R,i)] 
y/b x (t)b y (t)b z lt) 



A mt (R,0) 



(6) 



is contained in the scaling factors bj(t), scaled coordi- 
nates Rk = Rk/bk(t) and phase factor S(R(t),t) = 
S (t) + 2mJ2k[ R K t )/ b k(t)}db k (t)/dt. Now, we apply 
the time evolution results to an axially symmetric cigar 
shaped trap with uj x — u) y = lu±_ 3> u> z , and consider 
first the strongly interacting BEC limit. In the limit 
where £ = u) z /u>x <C 1, the scaling factor bj(t) can be ob- 
tained in powers of e, leading to b z (t) = l+e 2 [r tan^ 1 r — 
(l/2)ln(l + r 2 )] + 0(e 4 ) and b ± (t) = y/(l + r 2 ) + 0(e 2 ), 
where r = u±t is the dimensionless time. Therefore, ac- 
cording to Eq. © the aspect ratio of the cloud between 
the radial (L x , L y ) and axial (L z ) widths is given by 



L x {t) = b ± (t)u z = eb z {t) 
L z (t) b x (t)u x b±(t)' 



(7) 



In Fig [H the aspect ratio r xz = L x /L z is shown for 
e = 0.1 and e — 0.3, indicating that the anisotropy of 
the cloud is inverted during expansion and that the as- 
pect ratio reaches an asymptotic value significantly larger 
than one, similar to the expansion of Bose and 
s-wave Fermi [2l[ condensates. 

Notice that in the strong attraction (BEC) limit, 
fermions form tightly bound molecules and the internal 
degrees of freedom of fcrmion pairs do not play an impor- 
tant role. Thus, it is not surprising that the anisotropic 
nature of the p-wave interaction between fermions does 
not manifest itself in the behavior of the condensate. 
However, as one moves away from the BEC limit towards 
unitarity, the average fermion pair size increases, and the 
internal structure of pairs can dramatically change the 



condensate properties as the pair size becomes compa- 
rable to the inter-molecular spacing. In order to study 
this regime, we consider a similar case where a p-wave 
Fermi condensate is initially trapped in an axially sym- 
metric cigar-shaped potential with oj x — uj y 3> uj z . Fur- 
thermore, we impose a magnetic field applied along the 
x-direction (chosen as the quantization axis) to tunc 
through the Feshbach resonances. Since resonances in 
40 K are split [1, [H[ for different mi states, it may be 
possible to adjust the magnetic field such that fermions 
are paired in the mi = (p x ) state only. In this case, the 
p-wave interaction leads to the formation of p x symme- 
try pairs, which are more strongly correlated along the 
x-direction. As a consequence, the coefficients c* =0 in 
Eq. Q satisfy c x mi=0 > c y me=0 = c^ =0 , hence breaking 
the axial symmetry in the radial (x-y) plane. Thus, upon 
free expansion the cloud looses its axial symmetry which 
is initially imposed by the harmonic potential. 

In Fig. we plot the cloud aspect ratio in the radial 
plane r xy = L x /L y for the mi = state as a function of 
the p-wave scattering parameter l/{k' i F a p ) on the BEC 
side of the Feshbach resonance. In this plot, we choose 
the parameter ko/kp — 5, which is compatible with a 
dilute gas of 40 K with density n — 10 14 cm~ 3 . Notice 
that the anisotropy disappears in the BEC limit as it 
must, but becomes more evident with decreasing inter- 
action strength, leading to a 5% anisotropy near uni- 
tarity for sufficiently long times. Investigations on the 
BCS side [l/(kpa p ) < 0] and at unitarity require the 
inclusion of Landau damping which leads to the decay 
of Cooper pairs, as well as collisional dynamics of the 



normal flow [22j, |23j , which are beyond the scope of the 
present theory. In addition to the anisotropic expansion 
in the radial plane which is characteristic of the p-wave 
case, we also find an anisotropy inversion of the aspect 
ratio between the axial and radial directions, which is 
similar to the cases of Bose and s-wave Fermi conden- 
sates. 

The anisotropic expansion in the radial plane is bet- 
ter understood in the language of effective masses. In 
Eq. (HJ), the coefficient c l /d me can be absorbed in the 
effective mass Mj = md m Jd mi . In the BEC limit, 
Mi — > 2m is the mass of paired fermions. However, 
the effective mass becomes anisotropic for p-wave inter- 
actions as the system moves away from the BEC limit, 
such that for the p^-wave (mi = 0) discussed above, we 
have M x < M y = M z . Thus, it is easier to accelerate the 
cloud along the direction of lighter mass M x such that 
the cloud expands faster along the a;-direction than along 
the y-direction, hence breaking the axial symmetry. In 
contrast, if pairing occurs predominantly in the me = ±1 
channel (p y ± ip 2 -wave) , we have M x > M y = M z such 
that the cloud expands faster along the y direction, as 
shown in Fig. [2b- Thus, Figs. [2^ and[2jD help reveal the 
orbital symmetry of the p-wave order parameter. 

A clear distinction between the anisotropy inversion 
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FIG. 2: Aspect ratio in the radial plane r xy = L x /L y as 
a function of scattering parameter l/(k A F a p ) > at various 
time t (solid lines) for a p-wave Fermi condensate with (a) 
mg — and (b) mt = rtl orbital symmetry. Notice that r xy 
saturates to an asymptotic value (dashed lines) for r — > oo. 
Insets show the corresponding pair wave functions. 



sion of a cloud initially trapped in a harmonic poten- 
tial. We obtained the time evolution of the condensate 
on the BEC side of the Feshbach resonances for temper- 
atures close to the critical temperature. Starting from 
a cigar-shaped cloud with axial symmetry, we showed 
that the aspect ratio between axial and radial directions 
presents an anisotropy inversion during expansion, in- 
dicating a characteristic feature of condensate physics. 
Furthermore, we found an anisotropic expansion in the 
radial plane depending selectively on the magnetic quan- 
tum numbers mg = and mi = ±1 of the p-wave interac- 
tion. We also emphasized that this anisotropic expansion 
is a direct consequence of anisotropic p-wave interactions 
and occurs most clearly when the harmonic trap is com- 
pletely isotropic. In addition, we would like to stress that 
this anisotropic expansion should occur not only for p- 
wave, but also for any higher angular momenta (d-wavc, 
/-wave, etc..) as Feshbach resonances are approached 
from the BEC side. Finally, we proposed that the orbital 
symmetry of the order parameter for p-wave condensates 
can be directly probed through time-of-flight expansions 
of harmonically trapped clouds. Potential candidates for 
such experiments are 40 K, 6 Li, and 173 Yb. 

We would like to thank NSF for support (Grant N. 
DMR-0709584). 



induced by the harmonic potential and the anisotropic 
expansion induced by the p-wave interactions is found 
considering the time-of-flight expansion for a completely 
isotropic trap. In this case, for a p-wave Fermi conden- 
sate with either mg — of mg = H pairing, the cloud 
becomes axially symmetric during expansion and takes a 
cigar shape (for mg = 0) or pancake shape (for mg = ±1). 
This is in sharp contrast to the case of Bose and s-wave 
Fermi condensates, where an initially isotropic cloud re- 
mains alwyas isotropic throughout the expansion. 

Before concluding, we would like to stress that al- 
though Eq. (|4]) is obtained within the Ginzburg-Landau 
scheme for temperatures not too far below T c , our qual- 
itative results remain valid down to lower temperatures. 
In particular, the anisotropy inversion between axial and 
radial directions is just a direct consequence of super- 
fluid hydrodynamics, while the anisotropic expansion in 
the radial plane is due to anisotropic p-wave interactions. 
These two qualitative features are characteristic of p- 
wave Fermi condensates, and are not sensitive to tem- 
perature. Therefore, the qualitative aspects of our re- 
sults become more evident at lower temperatures, where 
a more sizable condensate fraction emerges (larger bi- 
modal distribution). 

In summary, we considered a single component p-wave 
Fermi condensate and studied the time-of-flight expan- 
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